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1 Introduction 

Active scalar problems are a wide class of research topics in fluid dynamics for which basic questions of 
existence and regularity pose challenging analysis problems that are often viewed as simpler ‘analogues’ of 
the famous Clay Math Prize Navier-Stokes problem. One subclass of such problems are the famous ‘vortex 
patches’ - exact L°° solutions of the incompressible inviscid two-dimensional fluid equations in which the 
scalar vorticity is the characteristic function of an evolving domain. The classical theory of general L°° 
weak solutions dates back to Yudovich ,56] in the early 1960s. However, the challenging problem of the 
long time regularity of the patch boundary was not settled until the early 1990s by Chemin m using 
methods from para-differential calculus. A geometric harmonic analysis proof was developed later by the 
first author and Constantin in [6]. 

More recent works have studied the dynamics of active scalars with a gradient flow structure. This 
opens up the possibility of using variational methods, including tools from optimal transport theory. The 
general problem has the structure 


d t p + div(pu) = 0, v = —VK*p (1.1) 

with many papers written to understand various subclasses of this problem El ii n m m m m nn 
on na nn [ 2 D m 123 nn m ho es eh eei gh uni m mi sa he E3 m es . This equation arises in many 
applications including materials science E3 ED * EU, cooperative control [ST] , granular flow [2.'j| ;25J [55], 
biological swarms [33] ill, [53] [51], vortex densities in superconductors [33 13 [Tj [32] S3 and chemotaxis 
[13 EH H3 QS]. Some of the recent literature has studied finite time singularities and local vs global 
well-posedness for both the inviscid case (11.11) [3 H3 HU 13 03 HD H3123 123 ED 123 [33] and the cases 
with various kinds of diffusion m na H3US3 uni Ea in multiple dimensions. The well-known Kellcr-Segel 
problem typically has a Newtonian potential and linear diffusion. For the non-diffusive problem ED, 
of particular interest is the transition from smooth solutions to weak and measure solutions with mass 
concentration. 

This paper concerns ED with the special case of K = N the Newtonian potential. This equation 
is exactly orthogonal to the classical 2-D inviscid incompressible fluid equations in vorticity form, in the 
sense that the velocity field is perpendicular to that given by the Bio-Savart law. Because of its gradient 
flow structure, the model makes sense in all space dimensions and we consider this problem in general 
dimension greater than one. We are interested in a special class of solutions called “aggregation patches”. 
These are particular L°° weak solutions in which the density p is a time dependent constant times the 
characteristic function of an evolving domain. Such solutions only exist for special kernels such as the 
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Newtonian potential. For this potential a previous paper [7] established a sharp C 7 and L°° regularity 
theory for solutions of (0> for K = N in general dimension. The proof generalizes the ideas of [J2] for 
C 7 solutions of the vorticity equation and [56] for the L°° theory. That work also developed a numerical 
method for computing aggregation patch boundaries and showed some interesting examples in both two 
and three dimensions. Of note is that these numerical solutions develop nontrivial geometric singularities 
at the blowup time - typically with mass concentrated along a “skeleton”-like structure. The simplest 
example with a trivial behavior is the collapsing sphere (or disk in 2-D) which due to symmetry collapses 
to a dirac mass at a single point at the blowup time. However elliptical initial data yield solutions that 
collapse onto a line segment and more complex initial conditions appear to collapse to a structure with 
branched arms. In an analogous fashion, the spreading solutions (backward time) were also studied both 
theoretically and numerically. In [7] the authors develop a rigorous theory proving L 1 convergence of 
the spreading patch solutions to an exact spherically symmetric similarity solution. However numerical 
simulations suggest that the weak L 1 convergence theory can not be made sharper due to the development 
of defects in the patch boundary in the approach to the similarity solution. 

These nontrivial dynamics open up the natural question of the regularity of the patch boundary for 
these aggregation patch solutions. In this paper we establish this result, working in Holder spaces as was 
done for the vortex patch boundary problem in fluids in E3E]. One key idea in [5j was a geometric lemma 
using cancellation properties of the gradient of the Biot-Savart kernel on half-disks and yielding a uniform 
estimate for the gradient of the velocity field with constants depending logarithmically on quantities 
measuring the smoothness of the boundary. We extend this logarithmic inequality to the multidimensional 
case involving even singular integrals and cancellation on half-spheres. It is worth mentioning that such 
kind of estimates have appeared before, without mention to the logarithmic dependence on constants, in 
connection with problems of classical analysis (see, for instance, m and US)- See also m , in which 
they appear in connection with the Muskat problem. 

A difficulty we have to confront in this paper that did not appear in the incompressibile fluids case 
is finding defining functions of the patch for non-zero times. The defining function for the initial patch 
transported by the flow is not smooth, because the field is not divergence free. We find a genuine smooth 
defining function adapted to our context, which leads to a commutator formula expressing the gradient 
of the velocity applied to the gradient of our defining function as a commutator of matrix valued singular 
integrals. A special Holder estimate in terms of the uniform norm of the gradient of the velocity field 
is derived. In addition, to prove our result we first develop the local existence and continuation theory 
for the patch boundary problem in general dimensions, which requires estimates of the transport map of 
the patch boundary in local coordinates. The technical apparatus needed at this point is more involved 
than the two dimensional case where one can parametrize by the circle, at least in the simply connected 
case. As a counterpart, we can work without additional effort on domains with holes and even on open 
sets made of finitely many pieces with disjoint closures. More details on the structure of the proof and 
information on the organization of the paper are given in the next subsections. 
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1.1 The main result 


Let d> 2 and let N(y) be the fundamental solution of the Laplace equation in R d . Thus N(y) = ^ log \y\ 
in dimension d = 2 and 


N(y) = -- 


1 


d > 3, 


(d-2)w d - 1 \y\ d ~ 21 

where ujd -i is the d— 1-dimensional surface measure of the unit sphere in ] 
equation 


1 . We consider the aggregation 


-£ + div(pu) = 0 
at 

v = —ViV * p 


( 1 . 2 ) 

(1.3) 


with initial data 

p(x, 0) = X Do 


(1.4) 


where Xd 0 is the indicator function of a bounded domain Dq C R d . We now fix 0 < 7 < 1 and take Dq 
to be a bounded C 1 " 1 " 7 domain (a domain with smooth boundary of class C 1+7 ; a formal definition will 
be presented in section HD- Then we have the following Theorem. 


Theorem 1.1. If Dq is a C 1+1 domain, then the initial value problem (11.21) . (11.31) and has a 

solution given by 

p(x,t) = 1 1 ^ XD t {x ), x £ R d , 0<t<l (1.5) 

where D t is a C 1+7 domain for all 0 < t < 1. 


As the proof shows, the preceding result also holds when Dq is a union of finitely many C 1+7 domains 
with disjoint closures. The conclusion is that D t is of the same type for all 0 < t < 1. It has been 
recently proved in [T] that the equation m flu has a unique solution in the weak sense for each 
initial condition Qo(x) in L°°(W 1 ) nL 1 (R d ). If the initial condition is the indicator function of a bounded 
domain Dq, then one has flUSD- In this case one speaks of aggregation patches, in analogy with the vortex 
patches for the vorticity equation associated with the planar Euler system (see [331 Chapter 8]). Thus our 
theorem solves the boundary regularity problem for aggregation patches. See 0 , [33] or |33] Chapter 8] 
for the analogous result for the vorticity equation in the plane. 

We describe a convenient reformulation of the problem that will be used throughout the rest of the 
paper. 

Set s = log(j3j), so that 0 < s < 00 if and only if 0 < t < 1. Define 


p(x, s) = (1 — t)p(x, t) and v{x, s) = (1 — t)v(x, t). 


Then, if the initial condition is dUD, to is equivalent to the transport equation 

dp(x, s) + \7 p[x, s) ■ v{x , s) = 0. (1.6) 

as 

The flows (or particle trajectories) in the time variables t and s are defined respectively by the ODE 


dX{x, t) 
dt 


v(X(x,t),t), X(x,0)=x 
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and 


—— = v(X(x, s), s), X(x,0) = x. 

as 

They are the same, in the sense that X(x, s) = X(x, t). Hence the solution of the transport equation (11.61) 
with initial condition p(x,0) = Xjj 0 (x) is 

p(x,t) = Xp^x), D s = X(D 0 ,s) = X(D 0 ,t) = D t = 

Dropping the tildes to simplify the writing and denoting again by t the new time s, 
problem (fra dra is equivalent to the non-linear transport equation 

dp(x,t) 

——— + Vp( x, t) ■ v(x , t) = 0 
v(x,t) = (—ViV * p)(x) 

with initial condition 

P(x, 0) = X Do (x). 

Theorem o can be then reformulated as follows. 

Theorem 1.2. If Dq is a C 1+7 domain then the initial value problem (11.71) . (11.81) and (ED has a solution 
given by 

p(x, t) = Xu t (a?), x G R d , tsR, 
where D t is a C 1+7 domain for all t £ R. 

The problem ED ED for d = 2 is similar to the vorticity equation for incompressible perfect 
fluids. The difference is that the velocity field in the vorticity equation is given by V^iV * p, which is an 
orthogonal gradient and, therefore, is divergence free. Instead the field (11.81) has divergence —p. 

1.2 Outline of the paper 

The proof of Theorem 1 1.21 is in two steps. First we look at the ODE giving the flow 

= v {X(a, t),t), X(a,0)=a, a € teK, (1.10) 

where the velocity field is 

v{x,t) = -(VN *XD t ){x), x€R d , D t =X(D 0 ,t). (1.11) 

Following Yudovich [56j (see j43j Chapter 8] for a modern exposition), the authors in [7] prove that 
(H.lOII - dl.llI) has a unique solution and that for each t the mapping a —> X(a,t) is a homeomorphism 
of onto itself satisfying a Holder condition of order /3(f), with /3(f) decreasing exponentially to 0 as 
t tends to oo. This does not use the smoothness of the boundary of the initial domain Dq and, in fact, 
it holds for an initial condition in L°°(R d ) D L 1 (R d ), with (11.111) modified appropriately. Assuming that 
D t is a C 1+7 domain for t in some time interval around t = 0, one can view equation (11.101) as an ODE 
in the Banach space C 1+7 (9Do,K d ). This ODE can be solved for short times by applying the Picard 
theorem in C 1+1 (8Dq, R d ) and thus one gets a flow of C' 1+7 -diffeomorphisms solving (11. 101) - (11 .111) . By 


we conclude that the 

(1.7) 

( 1 . 8 ) 

(1.9) 
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uniqueness of the Yudovich flow one concludes that the restriction of X(-,t) to 8Dq is of class C 1+7 on 
the surface 8Dq. In other words, the Yudovich flow is, for short times, of class C 1+7 in the directions 
tangential to 8Dq. This is discussed in section [2) Theorem 12.II provides the local existence result. One 
should remark that the statement of Theorem IQ includes a precise lower bound for the size of the time 
interval on which the solution exists, which will be used later on when dealing with long time existence. 
Showing that the Picard theorem can be applied is equivalent to various estimates, which are collected 
in Theorem 12.21 Its proof is presented in sections 0] and 01 One needs bounds for the action of principal 
value singular integral operators on Holder classes on smooth surfaces. The preliminary Section [3] gives 
two well-known lemmas in the precise forms that we need for the proof of Theorem 12.21 

The second step consists in proving that the Yudovich flow is of class C 1+7 in the directions tangential 
to 8Dq for all times. This requires a priori estimates for the quantities that determine the size of the 
local existence interval. The most relevant are those measuring the smoothness of the boundary of a 
C 1+7 domain. The C 1+7 smoothness of the boundary of a domain Dq is encoded in a defining function, 
that is, a function <f>o on class C 1+1 in R d , vanishing exactly on the boundary 8Dq and with non-zero 
gradient at each point of 8D 0 . By transporting <h 0 , that is, by setting ip(x, t) = <ho(-^ _ 1 (x, 1)) one obtains 
a function vanishing exactly on dD t , but with gradient \'(p(x,t) = Vd’o o VY r ~ 1 (x,t) which may have 
a jump at dD t , just as VX _ 1 (i,t). Thus tp(x,t) may not be a defining function for D t (and, indeed, 
it is not), contrarily to what happens in the case of the vorticity equation in the plane, for which the 
velocity field is divergence free. One of the difficulties that we have to overcome is finding a correct way 
of changing <f>o by means of the flow and still getting a genuine defining function <f>(x ,t) for D t . This is 
done in section [ 8 ] Once this is achieved, we need to get a priori estimates for the 7 -Holder semi-norm 
|| V<h(-, t )|| 7 on R d and for the infimum of |V<I>(x, t)| on dD t . The subtlest estimate is that of ||V<I>(-, i)|| 7 , 
which follows by bringing into the scene an appropriate commutator between a singular integral and a 
pointwise multiplication operator. This estimate is performed in section [3 Once the a priori estimates 
on the quantities determining the size of the local existence interval are available, the C 1 " 1 " 7 smoothness 
of dD t for all t, £ R follows readily. 

We close this section by showing that the transported defining function is already not smooth when 
d = 2 and the initial patch is the unit disc D$ = (x £ M 2 : |x| < 1}. The solution of the aggregation 
equation (mu-dill) with initial condition the characteristic function of the unit disc is 


POM) = Y~t XD (o>Vi=i)( x )’ xeR 2 , t< 1 . 


The field is 


and the inverse flow 


v(x,t) = -- 


1=7, |x| < y /1 -t, 

77 , |x| > y/l-t 


X i (x,<) = 




|x| < y/l-t, 
\x\ > y/l — t. 


y\ x \ 2 +t r- 

Take as defining function for Dq the function <^q(x) = |x| 2 — 1. Transporting ipo by the flow we obtain 

rh (M 2 - (!^ t)) , |x| < y/l-t, 


<p{x,t) = ip 0 (X (x,t)) = 


|x| 2 -(l -t), 


> y/T~t , 
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whose gradient has a jump at the circle |x| = yjl — t except for t = 0. To correct the jump one may take 

= (1 - t)xD t (x)ip{x,t) +XR2 \ Dt (x)<p(x,t), 


where D t = D( 0, y/1 — t). 


2 A Flow of C 1+1 Surfaces 


Given x £ with fixed d > 2 the cylinder with center x and radius r is 

C (x, r) = {x £ R d : |y' - x'| < r and |y rf - x d \ < r}, 

where we use the standard notation x' = (xi, ...,Xd- i) £ K d_1 . We say that D is a C ' 1+7 domain if for 
each x £ dD there exists r > 0 such that, after possibly a rotation around x, 


C{x, r) n dD = {y £ C(x, r) : y d = tp(y')} 


where p is a function of class C 1+1 in a ball B(x',r'), r' > r. In other words, the boundary of D is 
locally the graph of a C 1+7 function and thus a surface of class C 1+7 . A standard argument based on 
a partition of unity shows that if D is a C 1+1 domain then there exists a function $ £ (7 1 + 7 (R d ) such 
that D = {x £ R d : $(x) < 0}, dD = {x £ : $(x) = 0} and V$(x) ^ 0 for x £ dD. Such a function 

is called a defining function for D of class C 1+7 . Conversely, by the Implicit Function Theorem, if D has 
a defining function of class C 1+7 , then I? is a C 1+1 domain. There is a very useful quantity measuring 
the C ' 1+7 character of a domain, namely, 


q(D) 


l|V$ll7 

|V$| inf 


( 2 . 1 ) 


where 

|V$|i n f = inf{|V$(x)| : x £ dD} 


and, for each set E and each function / defined on E, we denote by ||/|| 7 ,b (or by ||/|| 7 , if there is no 
ambiguity on the domain of the function) the Holder 7 -seminorm 

ll/lk-E = sup{——j—- :x,y GE, x^y}. 

\x-yy 

The ODE providing the particle trajectories is (11.1 Of) - (1 1.11 fl . As we mentioned before, in [7] one proves 
that (11.101) has a unique solution and that for each t the mapping a —> X(a,t) is a homeomorphism of 
onto itself satisfying a Holder condition with exponent /3(t) decreasing exponentially to 0 as t tends 
to 00 . We call X(a,t) the Yudovich flow associated with the initial condition xd 0 - 

Assume that for t in an open interval containing 0 the restriction of X(-,t) to dDg is in C 1+7 (<9Do, R d ). 
Then dD t is a C 1+7 domain and we have 


v(x,t) = -(N * Vxr>J(x) = / N(x — y)n(y) da t (y) 

J dD t 

where n is the exterior unit normal vector to dD t and dat is the surface measure on dD t . If x = X(a, t) 
and we make the change of variables y = X(/3, t ) we get 


v(X(a,t),t) 



N{X{a , t) - X(J3, t)) (dX(/3, t)(T^ 3)) 


A •••A DX(p,t)(T d -im) MP), 
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where da is the surface measure on dDo, DX is the differential of X as a differentiable mapping from 
dD 0 into R d and Xi (/?),..., T<j_i(/3) is an orthonormal basis of the tangent space to dD 0 at the point 
/3 € dD 0 - The vector 

d -1 

DX(J3,t)(Tjm ( 2 - 2 ) 

3 =1 

is orthogonal to dD t at the point X(j3,t) and a different choice of the orthonormal basis Tj((3), 1 < 
j < d — 1, has the effect of introducing a ± sign in front of (EH). We may choose the Tj(fi) so 
that n(/3), Ti(/3),... ,Td-i(/3) gives the standard orientation of R d . Let D be the set of functions X e 
C 1+1 {8Dq ,R d ) such that there exists a constant /z > 1 for which 

\X(a) - X(j3)\ > -\a - /?|, a,PedD 0 . 

The smallest such /z is denoted by /z(X). Then X is bilipschitz and n(X) is the Lipschitz constant of the 
inverse mapping. It is clear that Q is an open subset of C 1+7 (cfDo, R d ). Given X € fl set 

. d-1 

F(X)(a) = / N(X(a) - X{0)) /\ DX^T^)) da(f3). (2.3) 

JdD 0 j=1 

Therefore X(a,t) satisfies the ODE 

= F (X(;t))(a), X(a, 0) = a, (2.4) 

which is called the “contour dynamics equation”. Our plan is to solve (12.41) for short times in the open 
subset of the Banach space C 1+7 (<9.Do,R d ). By uniqueness of the trajectories equation (ll.lOD - dl.llI) 
(see, e.g., [3j Theorem 3.7, p.128]), we conclude that the restriction of the Yudovich flow X(-,t) to dD 0 is 
of class C 1+1 for short times. In particular, dD t = X(dDo,t) is a surface of class C 1+1 for short times. 
In a second step we prove an a priori estimate which implies that dD t is of class C 1+7 for all times, thus 
proving Theorem 11.21 

Our estimates are most conveniently performed in terms of a particular norm defining the topology 
of C 1+1 (dDo, M d ). Let us momentarily drop the subindex ”0” and work with a bounded C 1+7 domain 
D. Since dD is a compact surface of class C 1+T there exists r = r(D) > 0 such that for each a € dD 
the set dD D C(a, r ) is the graph /3d = </?(/?') of a C 1+7 function t p (after a rotation around a if needed). 
The function 

X(p') = X(P, vtf)), P G B(a\ r) C R d ~\ 

is in C 1+7 (B(a' ,r)). Set 

v(X,a) = |Y(a)| + ||DX|| 00iB ( a / iT .) + ||DX|| 7iB ( a / ir ), 

where D is the ordinary differential of X and for a set E and a function / on E we denote by ||/||oo,.e 
the supremum norm of / on E. Finally set 

Ill'll 1+7 = Ill'll 1+7,0D = SU P v(X,a). 

a£dD 

Different choices of the local charts yield different but equivalent norms in C 1+7 (9D,R d ). 
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We discuss now two simple facts concerning the norm of C l+1 (dD, we just defined. The first 
is an estimate for the norm of the identity mapping I. In the local chart centered at a £ dD we have 
I(/3') = (/3 1 , ip(/3')), f}' £ B(a',r). Hence DI(/3') is a matrix with d rows and d—1 columns. The matrix 
formed with the first d—1 rows is the identity in R d_1 and the last row is (diip(p'), ..., d d -i(p(/3')) . Since 
we can assume that | V< J o(/3 / )| < 1, /?' £ B(a', r), by implicit differentiation we get (see the proof of lemma 
[HU] below) ||V<£>|| 7)jB ( a / )T .) < C d q(D). Then 

v(I, a) < |a| + C d + C d q(D). 

Let c be the center of mass of D and diam(H) its diameter. Then 

||^||i+ 7 , 9 d < diam(D) + |c| + C d + C d q(D). (2.5) 

The preceding inequality will be applied to dD t in dealing with long time existence. On the one hand, the 
center of mass is an invariant of the motion, because the kernel VX is odd. Without loss of generality we 
assume from now on that the center of mass of Do (and, consequently, of D t ) is the origin. On the other 
hand, we will obtain a priori estimates for diam(Z) t ) and q(D t ) (see (17.191) and (17.161) 1. We conclude that 
the estimate (ESP is good for our a priori estimates. 

The second fact we should discuss is how one estimates the Lipschitz constant of X £ C' 1+7 (9D,M d ) 
in terms of ||X||i_|_ 7 . Take points a,/3 £ dD. If \a — /3\ < r = r(D) we are in a local chart and then 
clearly |X(a) — X(/3)| < ||X||i+ 7 |a — /?|. Otherwise we estimate by the uniform norm and we obtain 
|X(a) - X(/3)| < 2||X|| 1+7 < |o - /9|. Hence 

\X(a)-X(/3)\ < (l + -^)\\X\\ 1+1 \a-/3\, a,/3 £dD. (2.6) 

Theorem 2.1. The initial value problem 

= F(X(.,t))(a), X{a, 0) = a, (2.7) 

has a unique solution X(a,t ) £ C' 1 ((—to, to), C 1+1 {dDo, R d )) and to depends only on d , q{Do),cr(dDo) 
and diam(Do)- 

This follows from the Picard Theorem for Banach spaces and 
Theorem 2.2. If X £ ft, then 

(a) 

F(X) £ C 1+ ^(dD 0 ,R d ) 

and 

||F(X)|| 1+7 < C 0 ^Xf d +\l + ||X||^+ 4 ), (2.9) 

where Co denotes a constant that depends only on d 1 q(Do) 1 u{Do) and diam(Do)- 

(b) X —> F(X) is locally Lipschitz on Cl: more precisely, 

||DF(X)|| <C 0 /x(X) 3d+8 (l + ||X||?^+ 7 ), Xefl, (2.10) 

where Co denotes a constant that depends only on d,q(Do),cr(Do) and diam(IDo). 
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Here by ||DF(X)|| we understand the norm of DF(X) as a linear mapping from C 1+1 (dDo, R d ) into 
itself. As the reader will realize, the precise form of the constant in the right hand side of (12.91) is crucial. 

We show now that Theorem 12.21 implies Theorem 12.11 The only point that requires further discussion 
is the size of the interval (— to, to) on which the solution exists. We need to find a ball B(I,p) C fl, of 
center the identity and radius p, so that F is Lipschitz in B(I,p) and we have an explicit bound for F 
on B(I,p). Lemma [iTTl gives for ro = t{Dq) the inequality rff 1 < 2q(Do). Take 

1 

P " i T i " x i x 
2(1 + 2 1+ 7< Z ( J D 0 )7) 

so that if X G B(I,p) then, by (12.61) . 

|X(a) - X(0)\ >\a-p\- \X(a) - a - (X(p) - /3)| 

>\a-P\ (l-(l+^)||^-/|| 1+7 ) 

\a — (3\ 

2 ’ 

that is, p(X) < 2, X € B(I,p). 

Clearly ||X||i +7 < ||/||i + 7 + p for X S B{I,p). By (12.101) F is Lipschitz in B(I,p) and by (12.91) 
||AT||i +7 < Cq 2 3d+2 (l + (||/||i +7 + p) 2d + 4 ) for X £ B(I,p). Therefore the solution of (12.71) exists in an 
interval (—to, to) with 

f > _ P _ 

°- Co 2 3d + 2 (1 + (ll^lli+ 7 + p) 2d+i ) ’ 

which is a quantity depending only on on d,q(Do),cr(Do) and diam(Z?o)- 


3 Two Lemmas 


To prove Theorem 12.21 we need two elementary lemmas on integral operators acting on Holder functions. 
These lemmas are well known but we prove them for the sake of the reader and because we need the 
precise dependence on various constants. 

Lemma 3.1. Let E be a measurable subset of R d_1 and assume that K : E x E —> R is a measurable 
function on E x E which satisfies 

\K(x,y)\ < (3.1) 

\K(x\,y) - K(x 2 ,y)\ < |a?i - x 2 \ ^ - y \d-i ' \^\ - xi\<\x\-y\/2. (3.2) 

Then for a constant C, which depends only on d and 7 , 


and 



K(x,y)f(y)dy 


00 , E 


< CAdiam(£) 7 ||/|| 00ii5 


(3.3) 



K(x,y)f(y)dy 


7 ,E 


<CA\\f\\^ E . 


(3.4) 
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Proof. Let D = diam(.E) be the diameter of E. We assume that diam(.E') < oo, otherwise (13.31) is trivially 
satisfied. If x £ E, then 


K{x,y)f{y)dy 




which is (13751) . For (E3) take two points Xi, X 2 € E and set 8 = 2\x± — X 2 \- Then 


(K(x!,y) - K(x 2 ,y)) f{y)dy 


< 


' Er\B(xi,8) 


\K{x!,y)\\f(y)\dy + 


f Er\B(xi,S) 


\K{x 2 ,y)\\f(y)\dy 


L 


\K(xi,y) - K(x 2 ,y)\ \f(y)\dy 

lE\B(x!,6) 

= Ii +12 + h- 

The first term is estimated as before : 

h < CAW/W^eST! = CA\\f\\ 00tE \ xi -x 2 \\ 

The term / 2 can be treated as I\ because E n B{xi,S ) C B{x 2 , §<5). For / 3 , by (13.21) . 


h < CAWfW^E \xi - x 2 


i 


dy 


(xi,s) ki -y\ d 7 


— < CAWfW^E \xx - a: 2 | 7 . 


□ 


Lemma 3.2. Let E be a measurable subset o 
E x E which satisfies 

\K(x,y)\ < 


-1 


and let K : E x E —»• R be a measurable function on 


— y\ d ~ l ’ 


\K(xi,y) - K(x 2 , y)\ < \xi - x 2 \ 


A 


\xi - y\ d ’ 


\xi - x 2 \ < \xi - y |/ 2 . 


(3.5) 

(3.6) 


Assume that b is a compactly supported bounded measurable function on E such that for a constant B 
one has 


’lloo.B + sup 

e>0 


K{x,y)b{y)dy 


<B, x £ E. 


I {y&E:\x-y\>e} 

Let f be a function on E satisfying a Holder condition of order 7 . Then 


K(x, y) (f(x) - f(y)) b(y) dy 


< CABdiam(Ey\\f\\ ltE 


and 


K(x, y) {f{x) - f{y)) b(y) dy 
E 

for some constant C depending only on d and 7 . 


<CAB\\f\\ ltE , 


7 ,E 


(3.7) 

(3.8) 

(3.9) 
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Proof. The inequality (13.811 is proved as in the previous lemma. Let us deal with (13.91) . Given x 3 ,x 2 £ E 
set S = 2\x\ — x 2 \. Then 


/ K{xi,y) (/(xi) - f{y))b{y)dy- / K{x 2 ,y) {f(x 2 ) - f{y))b{y) dy 
IE J E 

< [ \K{x 1 ,y)\\f(x 1 ) - f(y)\\b(y)\dy+ [ \K(x 2 ,y)\\f{x 2 )-f{y)\\b{y)\dy 

J EC\B(x\,8) J EnB(xi,5) 

[ K(x 1 ,y)(f{x 1 )~f(y))b(y)dy- f K(x 2 ,y)(f(x 2 )-f(y))b(y)dy 

IE\B(xi,S) JE 

= h+1 2 + h- 

The first term can be estimated readily by (1331) : 


h < CAB I 


,*[ 

J B(x i ,8) 


dy 


\x — 2 /| d_1_T 


——dy < CAB ||/|| 7 ,b \x\ - x 2 \ 


for some constant C depending only on d and 7 . For I 2 one only needs to observe that E fl B(xi,S) C 
E fl B{x 2 ,1 6). For I 3 we have, by (13.51) . (13.61) and (13.71) . 


I 3 < 


(f(xi) ~ / fa)) / K(x 1 ,y)b{y)dy 

J {y£E:\y-x 1 1 ><5} 


' {i/eB:|y-xi|>5} 


1 /(^ 2 ) - f(y )| \K(x!,y) - K[x 2 , y)\ \b(y)\dy 


< B\\f\\ Jt E \xi - Z 2 I 7 + CAB\\f\\ lt E \xi -X 2 \ f 
<CAB\\f\\ 1 , E \x 1 -x 2 \' y . 


dy 


{!/£R d - 1 :|x 2 -»|>|} I * 2 y\ d 7 


In the second inequality we applied (13.61) and then that \y — X\\ > ||y — x 2 \ and \y — x 2 \ > | for 
yGR d ~ 1 \B(x 1 ,6). □ 


4 Proof of Theorem 12.2L part (a) 


For convenience of notation we assume d > 3. The case d = 2 has a similar proof and can also be obtained, 
in the simply connected case, from the argument in [ 1431 Chapter 8 ] when the tangential derivative z a (a ', t) 
is replaced by the normal derivative iz a (a',t). 

Let X € fl and set ji = y(X). By (12.31) 

F(X)(a)= [ N(X(a)-X(P))G(P)da(P), a £ dD 0 , (4.1) 

JdD 0 


where G : 8Dq R d satisfies 

||G|| 7 ,aD 0 < C 0 \\X\\££ 

Since 3Dq is a compact C 1+1 surface there exists rg > 0 such that for each op £ 8Dq the part of 8Dq 
lying in the cylinder C(ao , 6 ro) is the graph ad = <p(a') of a function ip £ C 1+ ' y (B(a' 0 ,6ro)), after possibly 
a rotation around ag. We show in Lemma 16.41 below that we can take 


( 6 r 0 )- 7 = 2q(D 0 ) 


(4.2) 
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and that one has 


l|Vv3|| 7 ,S(ao,6r 0 ) < 2| ?(A))- 

Let if; G C£°(B(ao,3r 0 )) such that 0<^><1, ^ = lon B(ao,2r 0 )) and |VV’| < A/ro, A a numerical 
constant. Set F(X)(a) = Fi(a) + F 2 (a), a € dD 0 , with 

Fi(a) = [ N(X(a) - X{0)) G(/3) </>(/?) da(P), a € 8D 0 . (4.3) 

JdD 0 

Then F\ is the local part of the integral in ( 14 . 11 ) and F 2 the far away part. For a € dD$ D C(ao, 6ro) set, 
to simplify the writing, a = (ai,..., ad- 1 ) = a', so that (a, <p(a)) € dDo for a € B(ao, 6ro). Define 

Fji a ) = Fjfavia)), a£B(a 0 ,6r 0 ), j = 1,2. 

The function F 2 (a) is of class C°° in B(ao,2ro) and it is easily estimated in C 1+1 (B(ao,ro)) by taking 
gradient twice. The result is 

||*2||l +7)B (ao,ro) < 0) »(X) d (1 + ||X||f +7 ). 

The constant Cq in the preceding inequality contains explicitly the area of the surface a(dDo) and negative 
powers of ro, which can be estimated in terms of q(Do) by virtue of ( 14 . 21 ) . 

We turn now our attention to the more challenging term F\(a). To simplify notation set 

Z(a) = X(a, <p(a)), aeB(o 0 ,6r 0 ), (4.4) 

so that 

—^—r\a-b\<\Z(a)-Z(b)\<C 0 \\X\\ 1+1 \a-b\, a £ B(a 0 ,6r 0 ). (4.5) 

Define 

An estimate of the norm of Fi(a) in C 1+1 (B(ao, 3ro)) is equivalent to an estimate in this space of the 
function 

Tf(a) = J M(a,b)f(b)db 

where 

m = G(b, ipmw, <p(b)) (i + iv^i 2 ) 172 

is in G 7 (B(ao, 3ro)), has compact support in B(ao, 3t"o), and satisfies 

||/|| 7 ,B(a 0 ,3ro)<Go||^||^. (4.6) 

Passing to components we can assume that / takes real values. Our first task is to compute the distri¬ 
butional derivatives of T f. In view of the singularity of the kernel and the dimension of the space, which 
is d— 1 , we expect a singular integral of Calderon-Zygmund type to appear. That this is indeed the case 
is shown by the formula 

9jTf(a) = p.v. f -^-M(a,b)f(b)db, l<j<d-l, (4.7) 
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involving a principal value integral. The only difficulty in proving the above identity is to ascertain that 
the boundary term appearing in the integration by parts vanishes. If g £ C£°(B(a,ro)) then 


— [ Tf(a) djg(a) da = — lim [if M(a,b) djg(a) da) f(b) db. 
J e 'l° J \J\a-b\>s J ' 


Fix & and integrate by parts to get 
M (a, b) djg(a) da = i 


f JL 

\a—b\>£ d a j 


M(a,b) g(a) da + / g(a) M(a,b) rij(a) dcr(a), (4.8) 

J\a—b\=e 


J\a—b\>e 

where rij(a ) = (a, — bj)/\a — b\. To handle the boundary term in (14.81) . note first that since M(a,b) = 
0(\a-b\ 2 ~ d ) and a(dB(b,e)) = 0(e d ~ 2 ), we have 


lim / 

J\a—b\—e 


g(a)M(a , b)rij(a ) da (a) = g(b) lim 


e 4° J\a—b\=e 


M (a, b)rij (a) da(a). 


We now exploit the fact that rij is an odd function of £ = a — b to get 
M(a,b)n,j(a ) da(a) = 


0 


'\b—a\=e 


An elementary calculation gives 


f M(b + £,b)%da(t) 
l£l=c l?l 


J (M(b + £,&)- M(b - C, b )) da(0 


M(b + S,b)~ M(b - b)\ < ^p-yl Z(b + 0 + Z(b - 0 - 2Z(b)\ < 

which yields 

[ M{a 1 b)rij{a) da(a) = 0(e 7 ), 

J\a—b\=£ 

and, consequently, shows map . 

We now prove that the principal value operator in (14.71) maps boundedly C J (B(ao, 3ro)) into itself. 
The strategy is as follows : if there were the derivative with respect to bj in the kernel of the operator in 
dUD we could try an integration by parts. We show that a sort of commutator changing the derivative with 
respect to aj into one with respect to bj has a kernel with extra smoothness and thus the corresponding 
commutator operator satisfies the C 7 estimate we are looking for. Set 


C{a,b) = (2 -d)- 


(Z(a)-Z(b)).(J-Z(a)-J-Z(b)) 


\Z(a) — Z(b)\ d 


We then have 


d „„ „ , n JZ(a)-Z(b))-£Z(a) 
daj M M ( 2 d ) \Z(a)-Z(b)\ d 


= C{a , b) + (2 - d) 


( Z[a) - Z(b)) • ZLZ(b) 
\Z{a) — Z(b)\ d 


= C{a,b) - —M(a,b) 

dbj 
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To show that the operator 

Cf(a) = J C(a,b)f(b)db (4.9) 

is bounded on C 7 (B(ao, 3ro)) we appeal to lemma [3721 Remark that 

f)7 r)7 

C(a,b) = K(a,b)-(—(a)-—(b)) 

with K(a,b) a kernel which satisfies the hypothesis (13.11) and (13.21) of lemma 13.21 with constant A = 
Co\\X\\\ +1 n(X) d+2 . To apply lemma [3~51 we need to check that 


f\b—a\>e 

For that write, for a £ B(ao , 3ro), 


K(a, b) f(b) db 


< C, e > 0, a£B(ao,3ro). 


/ K(a,b)f(b)db = K(a, b)(f(b) — f(a)) db + f(a) 

J e<\b—a\ J e<\b— a\<6ro 

= 1(a) + 11(a) 

The term I (a) is estimated straightforwardly by 


£<|b— a|<6r 0 


K (a, b) db 


\I(a)\<C 0 \\X\\ d 1+ ^(X) a 


J\a—b\<&r 0 I® 

For 11(a) we use a “pseudo-oddness” property of the kernel. We have 
f Z(b) - Z(a) 


-—-db = C 0 \\X\\{ +ll ji(Xy 


£<|a—b|<6r 0 \Z(b) — Z (a)\ d 

Z(a + q)-Z(a) 
:|4|<6r 0 \Z(a + 0 — Z(a)\ d 


db 


dH 


Z(a + £)-Z(a) + Z(a-£)-Z(a) 


2 Je< |£|<6ro VI Z(a + £) - Z(a)\ d I Z(a - 0 - Z(a) 


d( 


[ (Z(a + £)-Z(a)) ( 

^ 2e<|{|<6r 0 \ 


\Z(a + 0~Z(a)\ d \Z(a-(()-Z(a)\ d 


dti 


L 


Z(a + a) + Z(a-0-2Z(a) 


|?|<6r 0 I Z(a~a) — Z(c 


da 


The elementary inequality 


\z\ d -\w\ d 


<d sup (\z\ d 1 i ,\w\-’)\z±w\ 
0<j<d-l 


provides an estimate for the first term above and the second is estimated straightforwardly. We finally 
obtain 

\H(a)\ < Co Ll(X) 2d (1 + ||X||*£) ||/||oo,S(ao,3ro) 

<Con(X) 2d (l + \\X\\l d + y. 
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The constant of the kernel I\(a,b), as in (13.11) and (13.21) . is less than C 0 fi(X) d+2 ||X|| 3 +7 . Therefore 
lemma 13.21 yields 

||C'/|| 7> B (a0l 3r 0 ) < Co /i(X) M+2 (1 + ||X|0. (4.10) 

It remains to estimate the operator 

Uf(a) = p. v. j -^-M(a,b)f(b) db 

on C 7 (B(ao, 3ro)), for 1 < j <d — 1. Take a function \ 6 C“(I3(ao,4ro)) such that 0 < y < 1, y = 1 on 
B(a 0 ,3r 0 ) and |Vx| < C 0 /r 0 . Then / = /%. We have 

Uf(a) = j A. M(a,b)(m-f(a))x(b)db 

+ f(a) p,./ M(a,b)x(b)db 
= U 1 f(a)+f(a)U 2 (a). 

Integrating by parts and noticing that, as before, the boundary term vanishes we get 

/ d 

M{a,b)—x(b)db 1 a € H(a 0 , 3r 0 ). 

Thus, by (031) and ||'VxlU,B(a 0 ,3r 0 ) < CoTq 1 , 

\U 2 (a)\<Con(X) d ~ 2 , a6B(fl 0 ,3r 0 ), (4.11) 


and by lemma l3Jl 

II C 2 11 7 , B(ao,3ro) < C 0 fJ-(X) d 1 1111 1 +0 ,- 

Here the constant of the kernel M(a, b) has been estimated by Co /x(X) d_1 ||X|| 2 +7 . By (14.61) 

ll/^2|| 7 ,B(ao,3n,) < Co n{X) d ~ l WXWi+i. (4.12) 

For U\f we apply lemma [3~2l The kernel of the operator U\ is d/dbjM(a,b), whose constant turns 
out to be not greater than Co fJ,(X) d+2 (1 + ||X||f +7 ). Tacking into account (14.61) and (14.111) lemma HOI 
yields 

l|Ci/|| 7iB(ao , 3ro ) <C 0 ^X) 2d (l + \\X\\ d ^). (4.13) 

Combining (14.121) and (14.131) 

l|C/|| 7l B(a 0 ,3r 0 ) < C 0 »(X) 2d (1 + ||X||*£). (4.14) 

By (14.141) and (14.101) we finally obtain 

II W(a)|| 7 , B(o0i3ro) < Co M (X) 3d+2 (1 + ||X|| 2 t+ 4 ), 
which completes the proof of (EH). 
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5 Proof of Theorem 12.21 part (b) 

It is enough to show that for X £ Q and H £ C 1+7 (<9.Do, R d ) 

DF(X)(H) = 4tF(X + XH) 

(A A X—0 

satisfies 

\\DF(X)(H)\\ 1+1 < C 0 n(X) 3d+s (1 + ||X||“+ 7 ) ||H|| 1+7 . (5.1) 

To prove this we may assume that ||K1 ||i + 7 = 1. We first compute 


DF{X)(H) = -ttF(X + X H) 

uA 


A —0 


d_ 

dX 


« d— 1 

/ N((X + XH){a)-(X + XH)(pj) /\ D(X + XH){T 3 {fd)) da(/3) 
A=0 JdD 0 V ' j=1 


d—1 


i a d 0 


N(X(a) - X(J3)) £(-l y-'DHiTjm /\ DX{T k {t3)) da((3) 


i=i 


4-1 


/ VX(X(a) - X(P)) ■ ( H(a ) - (0)) A D A)(T 3 m MP) 

ldD 0 j= \ 


= A(a ) + B(a) 

Consider first the term 4(a). This is a sum of d— 1 terms, each of which looks like the function F(X)(a) 
in El. The only difference is that in A(a) one of the factors DX(T 3 )(/3 ) has been replaced by a vector 
of the type DH(T 3 )(/3). Then the estimate of A(a) in C 1+1 (DDq, R d ) is performed in exactly the same 
way as we did in the previous section for F(X). There is only one difference, namely, that in the bounding 
terms one of the factors ||X||i_j_ 7 should be replaced by ||ff||i +7 = 1. Thus 

||^|| 1+7 < C'oM(^) 3d+2 (1 + lbY||?^.+ 3 ). 

The term B(a) is slightly different because of the presence of the factor H(a) — H(/3) in the kernel, 
which compensates the higher singularity of VX(X(a) —X(/3)). The structure of the argument is, however, 
the same. One performs the splitting into local and far away parts, as in (Id.3D . The local part, which is 
the most difficult, can be written in local coordinates a = (ai, ...,a,d-i) as 

Tf(a) = jM(a,b)f(b) db 

where / is a scalar function satisfying the estimate (14. 6D . and the kernel M(a, b) is given by 


M(a,b) = XN(Z{a) - Z{b))-{h{a) - h(b)), 

Z(a) = X(a,ip(a)) and h{a) = H(a,<p(a)). The function Z satisfies (14.51) and ||h|| 1+7iB ( a0i3r , 0 ) < Co- As 
before, the boundary term vanishes and we have 

r f) 

djTf(a ) = p.v. / ——M(o, b)f{b) db, l<j<d-l. 
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We express this operator as a commutator minus an operator with kernel dj/dbjM(a,b). Recall that the 
commutator gives the worst constants. The kernel of the commutator is 


C(a,b) = V 2 1V (Z(a) - Z{b))(d j Z{a) - djZ(b))-(h(a ) - h{b)) + X7N(Z{a) - Z(b))-{djh(a) - d 3 h{b )) 

= K{a , b)(djZ(a) - d 3 Z(b)) + X7N(Z{a) - Z(b))-{d 3 h(a ) - d 3 h{b)), 

where the second identity defines the matrix K(a , b). The operator given by the kernel in the second term 

WN{Z(a) - Z(b))-(djh(a) - d 3 h{b )) 

is estimated as we did in the previous section for flUl). The worst constants appear in estimating the 
operator with the kernel 

K{a,b){d j Z{a) - d 3 Z{b)). 

We follow closely the argument for the estimate of GH). The step that gives the largest constant happens 
when dealing with the quantity 

/(a) f K(a,b)db, a£B(a 0 ,3r 0 ). (5.2) 

J e<\b—a\<6ro 

The pseudo-oddness property of K gives 

\[ K(a,b)db\<C 0 \\h\\ 1+1 v(X) 2d+4 (l + \\X\\i+ 4 1 ), aeB(a 0l 3ro), 

J e<16—a|<6ro 

which, combined with GH), yields the upper bound 

C 0 /x(X) 2d+4 (l + ||X||^+ 3 ) 

for the norm of the matrix in & The constant of the kernel K(a,b) in applying lemma f- 1.21 is 
Cq n(X) d+i ||X|| 4 +r Thus lemma [372] finally gives 

||T/|| 1+7 <C 0 /i(X) 3d + 8 (l + ||X||^+ 7 ) 

and 

||Z)F(X)(£0|| 1+7 < Co MX) 3d+8 (1 + imi“+ 7 ), 
which is (15.Ill , because we are assuming that ||if||i +7 = 1. 

6 The logarithmic inequality for || Vv(*, t) ||oo 

Fix a bounded C 1+7 domain D and write 

v(x) = — V1V * Xd(x). (6.1) 

In section we will establish a logarithmic estimate for UVuHoo that will be needed to get long time 
solutions of the problem (fTTTl) . (OH and USD- 
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Lemma 6.1. Let x dD, and let e = e[x) = dist(;r, SD). Then for 1 < j,k < d, the vector v = 
(v 1 , v 2 ,..., v d ) satisfies 


dyj < , 

0*0 = 

OXk 


Xj Xfo 


Vd -1 F 


d +2 


and 


dv\ , l v , , /I M 2 -da;? 

^-W = - 5 x D W-l— —j^p+5— 


0*0, j + k 

(6.2) 

* %D\B(x,e) 0*0 

(6.3) 


where uid —1 = <r(§ d 1 ) and the derivatives are distributional derivatives. 
Proof. Suppose j ^ k. By (16.11) and Green’s theorem, 

(xj - Vj)(xk - yk) 


— (x) = d lim / 

dx k y u d -i ^°J Dr Uv_x\>n\ \x-y\ d + 2 


dy, 


but for 0 < r\ < e{x) 


I Dn{\y-x\>ri} 

(xj - yj)(x k - yk) 


dri<\y—x\<e \ x V I 

That established & and the proof of (16.31) is similar. 


d+2 


dy = 0. 


□ 


Notice that the principle value kernels in m and (lOll have the form 

= x ^°> ( 6 - 4 ) 

where 

(i) is homogeneous of degree 0, SI (a;) = S2(A), 

(ii) SI is even, il(—x) = Sl(x), 

(hi) SI G £7 1 (R ,, \{0}) ! 

and 

(iv) J gd _ 1 fl(x) da(x) = 0. 

As we mentioned before, a C 1+1 domain D has a defining function, that is, a C 1+7 function $: —> R 

such that and £) = {$< 0}, dD = {<h = 0} and V$(:r) 0, x G dD. We set 


IV<E>| in f = inf |V$(x)| 

x£oD 


and 


||V$|| 7 = sup 


|V$(X!) - V$(s 2 )| 


l 'y | | 

Xl^X2G^- d I**! x 2y 

Theorem 6.2. Let K satisfy (16.11) and (i)-(iv) and let D be a C 1+1 domain with defining function $. 
Then 


K(x — y)X D (y) dy 


< 


^i + iog+Ani^lMk^, I£ l d , 0 < e, (6.5) 
l\y-x\>e 7 V V |V$|inf// 

where Cd and C(S7) are constants that depend only on d and S7 respectively, and log + x = maxjlogx, 0} 
is the positive part of the logarithm. 
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Corollary 6.3. Ifv(x) is defined by (16.11) and & is a defining function for the bounded C 1+7 domain D, 
then 


C 

|Vw||oo < — 

7 


1 + log 4 



\m\y X\ 
|V$|in fJJ ■ 


( 6 . 6 ) 


The Corollary is immediate from (16751) . (El and El- 
We need a lemma. 


Lemma 6.4. Let D be a C 1+1 domain with defining function $. If S > 0 satisfies 

j ,[^SJ 2< i 

|V®U “ 2’ 


then for each x £ dD, dD U B(x,5) is, after a rotation around x , the graph of a C 1+7 function and 
D D B{x, S) is the part of B(x, S) lying below the graph. 


Proof. Assume, without loss of generality, that x = 0 and that V$(0) = (0,..., 0, ^^(O)), ^$(0) > 0. 
Take two points p,q £ dD U B( 0, 6) and set p = (x 1 , Xd) with x' = (a;i,..., Xd- i), and q = (y’, yd) with 
y' = (yi, ■ ■ •, 2/d—i)- Then 


0 = $(p) = $(0) + V4>(0) • p + E(p) = |V$(0)|x d + E(p) 


and similarly for q. Subtracting and taking absolute value 

|V$(0)||* d - y d | = | E(p) - E(q)\ < ||V$|| 7 5 7 (\x' - y'\ + \x d - y d \) 

and thus 


I x d - Vd\ < ||V4> 
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which yields 


\x'-y'\ ~ |V$| inf 
| x d - y d 


<5 7 1 + 


\x d - Vd\ 


\x' — y’ 1 ' 


< 1. 


\x’ - y '| 


This says that dDL)B(x , S) is the graph of a Lipschitz function x d = <p{x'), with domain an open subset U 
of {x' £ R d_1 : \x'\ < <5}, satisfying iV^x 7 )! < 1, x' £ U. By the implicit function theorem ip is of class 
C 1+7 on its domain and this completes the proof of the Lemma. 

Notice that U contains the ball {x' £ R d_1 : \x'\ < <5/2 1 / 2 }. We need also the estimate 


By implicit differentiation 


IV^COI < (2rf) 1/2 r 7 , \x'\ <r< W. 


x',p{x')) 


(6.7) 


and so 


which gives 


\dj<&(x',p(x'))\ < \d d ^(x',<p(x'))\, l<j<d-l, 


|V$(a/, tp(x'))\ < d 1/2 \dd$(x',p(x'))\. 


Since 9j$(0) =0, 1 < j < d— 1, 


|V4>|inf 


which completes the proof of El- 


□ 
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Proof of Theorem \6.2\ Assume first that x £ dD. Take S > 0 such that S ' 1 = L and set r] = d/2 1 / 2 . 
Let e satisfy 0 < e < r/. Then 


[ K(x - y)X D (y)dy = [ K(x-y)dy+ f 

J\y—x\>e J Dn{e<\y—x\<n\ JL 


Dn{e<\y—x\<ri} 


K(x - y) dy 


Thus 


\h\ < f 

J L 


and 


f\ D \ 1/d dr 

|-^21 — ||^||oo^d—1 / — 11^| |oo^d—1 log 

Jr] ^ 


Dn{\y-x\>\D\i/ d } \ x — V I 

D\ X ' d 


/Dn{i)<|x-i/|<|r>| 1 /' i } 

+ K(x-y)dy = h+h +h. 

J Dr{\y—x\>\D\ 1 / d } 

l ^^dy <%^| = |MU 


\D\ 


y 


)<ftC(!!) (l + log+^rf^t)). 


If rj > \D\ 1 / d , then we let I 2 = 0 and this brings in again the positive part of the logarithm. 

Let us turn to I\. Assume that x = 0 and that we are in the situation discussed in the proof of 
Lemma 16.41 Taking polar coordinates we get 


\h\< 


[ [ 17 (w) da(ui) — 
J e J A(r) ^ 


where A(r) = {cj G 5 d_1 : rco G D}. Let H stand for the half-space {x G : Xd < 0}. Since is even 
and has zero integral on S' d_1 , we conclude that the integral of 17 on the hemisphere S^ -1 D H is also 
zero. Hence 

/ 17 (uj)da(uj)= / 17(w) da(uj) — / 17 (uj)da(uj), 

JA(r) JB(r) JC{r) 

where B(r) = {w £ S' d_1 : roj £ D\H} and C(r) = {w £ S' d_1 : roj £ H\D}. Let us proceed to estimate 
the integral on B{r) (the integral on C(r) is estimated similarly). 

For some absolute constant Cg (which can be taken to be n/2) one has, by (16.71) . 

a{B(r)) < C 0 sup \<p(x')\- < C 0 sup |V^(x')l < C 0 (2d) 1/2 r 1 . 

r \x'\<r |V$| inf 


\x'\<r 


Therefore 
rv 


dr 


/ 17(w) da(uj) 

e J B(r) r 


< INloo f\(B(r))- < ||17|| 00 C' 0 (2d) 1 / 2 M^i rf < C d C(17) 1 

J e ' 


|V$| in f 7 


7 


which completes the proof for 0 < e < rj. If y < e, then I\ = 0 and I 2 and I 3 are estimated as before. 
Let us assume now that x dD. Let eg denote the distance from x to dD. If e < eg, then 


f K{x - y)X D {y)dy = f K{x-y)X D {y) 

J\v—x\>e J\v—x\>en 


'\y-x\>e J\y-x\>e 0 

and so we can assume that eg < e. Take xq £ dD with \x — xq\ = eg and define 


dy 


A = 


[ K(x- y)X D (y) dy - f K ( x 0 - y)X D (y) dy. 

J\y-x\>e J\y-x 0 \>2e 
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Then A = Ai + A 2 , where 


and 


We then have 


Ai = [ (k(x- y) - K(x 0 -y)\x D (y)dy 

d\y-x 0 \>2e^ ' 


A 2 = / K(x - y)X D (y)dy. 

J B (xq ,2e)\B (x ,e) 


|a 2 [ < Halloo |jg(x °; 2f)l =Cdll»||c 


and, by a gradient estimate, 

|Ai | < C(f2) \x - x 0 | [ 

J 2, 

which completes the proof of Theorem 16.21 


dy 


2e<\y—xo\ \ X ^ol^ 1 


<c d c(n), 


□ 


7 Global Existence 

We prove in this section that the Yudovich flow X(a,t ) solving the ODE (11.101) and (11.111) is smooth 
in the directions tangential to 8 Dq for all t £ K. More precisely, the restriction of X(-,f) to DDq is in 
C 1+1 (dDo, R d ) for all times f 6 R. In particular, dD t is a domain of class C' 1+7 for t £ R. The local 
existence Theorem 12.11 shows that X(-,t) is in C' 1+7 (9Do,K d ) for t £ (—T, T), where T is the small time 
given by the Picard Theorem. Assume that T is maximal with the property that the solution X(-,f) is 
defined in (—T, T). Our goal is to prove the a priori estimates on (—T,T) which will let us to conclude 
that indeed T = 00 . 

It is enough to prove that dD t is a domain of class C 1+7 for all f G R. In fact, if this is true, then we 
have 

||A||i +7i aD 0 < 00 , teR. (7.1) 

Otherwise, let T be a maximal time so that (17.11) holds for t £ (—T,T). Taking Dt or as initial 
domain at time T or —T (not at time 0 !) in Theorem 12.II we contradict the maximality of T. 

To show that dD t is a domain of class C 1+1 for all t £ M take a defining function <f>o for Dq. Then 
$0 is in C 1+7 (K £i ), D 0 = {<E> 0 < 0}, dD 0 = {<t> 0 = 0} and Vd>o(a:) /0,i£ dD 0 . Consider the equation 
d 

— $(x,t) + V$(x,t)-v(x,t) = — $(x,t)xD t (x), x £ R d , t£R , (7.2) 

with initial condition $(x, 0) = $o(a ; )- Then 

<h(x,f) = $ 0 (A _1 (a;,t)), x£R d \D t , 

and 

$(x,t) = e~ l $ 0 (A _1 (a;,t)), x £ D t , 

where, for a fixed time t, X~ 1 (x,t) is the inverse of the mapping x —> X(x,t). Notice that &(x,t) is 
continuously differentiable in the open sets D t and \ D t , but that V<l>(a;, t) may, a priori, have a jump 
at the boundary of D tl just as VT _1 (i,t) or X7X(x,t). We will show in the next section that $(x,t) is 
of class C 1+7 (R d ) and thus a defining function for D t . We use this fact freely in this section. 

The a priori estimates we need are collected in the following statement. 
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Theorem 7.1. Let 4>(-,f) t/ie defining function for D t determined by (17.21) . Then 


and 


l|V$(-,t)||oo < ||V$(-,0)|| oo exp f (1 + ||Vv(-,s)||oo)rfs, 

Jo 

|V$M)| inf > |V$(-,0)| inf exp(- [ (1 + IIVy(-,s)||oo) ds) 

Jo 

||V4>(.,f)|| 7 <||V4>(.,0)|| 7 exp(C /V+ l|Vw(- ! s)|| 00 )rfs). 

Jo 


Proof. Taking derivatives in m we obtain that the material derivative of V$ is 

-g(V$) = -Vr.(V$)^ XDt V<I>. 

We have used here that, since $(x, f) vanishes on dD t , 

V\D t (x) = $(x,t) n{x) da(x) = 0. 

By CU) 

|V$(*,i)|<|V$(®,0)| + [ (l + ||Vu(-, s)^) ||V$(-, s)^ ds 

J O 

and US follows from Gronwall. 

For (17.41) take x € dD t . Then 


21 log |v«(i, t)l = |v ^(L,t)p 




> -- 


1 


|V$(s,i)| 


D 

~Dt 


(V$(x,f)) 


>-(i + IK,f)|| 00 ) 


and so 


|V$(x,i)| > |V$(x,0)| exp J o (1 + ||Vu(-,s)||oo)ds^ , 


which yields m at once. 

For (17.511 we need two lemmas. The first one is an elementary remark. 

Lemma 7.2. If D is a bounded (7 1+7 domain and & is a defining function for D, then 

Xd(x)<&(x) = VN * (xuV<h)(x), x € R d . 


(7.3) 

(7.4) 

(7.5) 

(7.6) 


(7.7) 


Proof. On one hand, the functions in either side of (TTH) are continuous functions. On the other hand, 
we have the distributional identities 

ViV * (xdVcE*) = VAT * V(xd$) = AN * xd& = Xd&- 


□ 
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Denote by HN the distributional Hessian matrix of N. If i ^ j the entry in HN corresponding to the 
*—th row and j— th column is the distribution 


p.v. dfjN(x) = p.v.-- 


XiXj 


V d -1 |x| d + 2 ’ 

while the diagonal term corresponding to the indexes i = j is the distribution 

„ 2 »r, x 1 r 1 \x \ 2 — dxf 1 - 

p. V. d u N(x) + - d 6o = p. v. — -™r- +-So, 


(7.8) 


(7.9) 


where <5o stands for the Dirac delta at the origin. In (17.811 and (17.911 the second order partial derivatives 
of N are taken pointwise for x ^ 0. Hence 


HN = p.v. X/ 2 N + -Io, 


(7.10) 


where Iq stands for the diagonal matrix with Sq in the diagonal. In the next lemma we establish a 
commutator formula which is crucial in what follows. In the statement below N 2 N(x), x ^ 0, stands for 
the dx d square matrix with entries the pointwise partial derivatives dijN(x). The integrand in the right 
hand side is absolutely integrable because the vector V$ satisfies a Holder condition of order 7 . 


Lemma 7.3. 


^-(\7<S>{x,t)) = [ V 2 N(x-y)(N$(x) - V$(y)) dy, 

m JD t 


x e 


Proof. Since v = —NN * \D t , taking gradient we get 


1 


Vw = -HN * X D t =- p. v. V N * XD t - ~,XD t I, 


with I the identity matrix. Identities (17.61) and (17.121) yield 


-^-(V$) = (p. v. V 2 iV * XD t ) (V$) + ^Xd, VT - X Dt v$. 


Taking gradient in (ED 

XD,V'h = HIV* (xD t V$) = p.v. V 2 1V* (X£» t v$) + -XD t V$. 

d 

Combining (17.131) and (17.141) completes the proof of the lemma. 


(7.11) 

(7.12) 

(7.13) 

(7.14) 

□ 


Lemma 17.31 yields the a priori estimate ED exactly as in [ 6 ] or [43] . Theorem 17.II is then proved. □ 

Inserting in (16.61) (the logarithmic estimate of || Vu(-, t)||oo in terms of q(D t )) the a priori estimates 
ED and E we get 

||Vv(a;,£)||oo < C + C [ ||Vv(-,r)||oo dr, 

Jo 

which yields by GronwalPs inequality 


||Vn(a;,t)||oo <Ce 


ct 


—T <t <T. 


(7.15) 
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Upon introducing this in (17.31) and (17.dl) we conclude that we have the double exponential estimate 


q(D t ) = 


l|Vd>(-,t )|| 7 


< C exp(e Ct ), —T < t <T. 


|V$(', t)\i n f 

For fixed t, the flux a —> X(a, t) is a bilipschitz homeomorphism. In particular, we have 


||V.Y(a,t )|| 00 < exp (^J ||Vw(-, s)||oo ds 


and so 


a(dD t ) = f I A DX{P,t){T 3 m dcr([3) 

JdD 0 1 j=1 

< [ {d~l)^\\DX{;t)\\Vd^) 

JdD 0 

< (d- 1)5 exp (^{d- 1 ) J ||Vu(-,s)||oods^ a(dD 0 ). 


Hence, by (17.151) . 


a(dD t ) < (d— l ) 2 cr(dD 0 ) exp (C e Ct ), —T <t <T. 
A similar estimate of the diameter of D t follows from (17.151) and (17.171) : 

diam(H t ) < diam(H 0 ) exp (C e Gt ), —T <t <T. 


(7.16) 
(4.47), p.149] 

(7.17) 


(7.18) 


(7.19) 


We can combine the estimates (17.161) . (17.181) and (17.191) with Theorem 12.11 to show that D t is a C 1+1 
domain for all t & R and thus complete the proof of our main result. For that, assume that T < oo is the 
maximal time for which D t is a C 1+7 domain for all t £ (—T,T). Since the size of the interval in which 
the local solution of Theorem rrn exists depends only on the quantities q(Do) 1 a{dDo) and diam(Zlo) 
and these quantities are uniformly bounded for D t with t £ (—T, T), we can apply Theorem 12.11 with 
initial condition D Tq at time To (not at time 0!) with To sufficiently close to T so that the new interval of 
existence goes beyond T. This contradicts the maximality of T. 

Notice that we have concluded the proof without proving any a priori estimate for \\DX(-, t)|| 7 ,a_D 0 - 
This is the reason why we took particular care in estimating the local time of existence in terms of 
quantities related to the smoothness of the initial domain D$. A final remark on inequality (17.191) is in 
order. One can easily prove the much better estimate diam(D t ) < Cd diam(Do), with Cd a dimensional 
constant, using the straightforward potential theoretic estimate \v(x,t)\ < Cd\D t \ 1 ' d , x £ t £ R, 
and the fact that \D t \ = e~ t \Do\. 


8 The gradient of $ has no jump 

In this section we prove that the function $(a :,t) defined by (17.21) is of class C 1 + 7 (R d ). 

Let X(a,t) be the Yudovich flow (11.101) and (11.111) . The initial domain Dq is bounded and has 
boundary of class C 1+7 , 0 < 7 < 1. By Theorem 12.II we know that for some T > 0 


X(a,t) £ C 1 ((—T, T), C 1+1 (dD 0 , R d )). 


( 8 . 1 ) 
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Fix a time t, 0 < t < T. By (16.61) || Vu(-, t)||oo is finite and controlled by the constant q(D t ) describing 
the C 1+1 character of dD t . In view of (18.11) this constants are uniformly bounded for |r| < t. Hence, by 

(EU), 

||Vu(-,r )|| 00 < q, |r| < t, (8.2) 

for a positive constant q depending only on t. We know from Lemma 16.11 that the entries of the ma¬ 
trix Vv(-,t) are given by singular integrals with even kernels plus a scalar multiple of X o t ■ By the main 
lemma of 32] Vu(-, t) satisfies a Holder condition of order 7 on each of the open sets D t and R. d \ D t , in 
spite of having a jump at dD t . Again the constants of this Holder conditions are controlled by the C 1+7 
character of dD t . Therefore, for some other constant q depending only on t, 


||Vu(a;i,T) - Vu(x 2 ,r)|| < q |xi - x 2 \ J , |r| < t, 


(8.3) 


provided xi,x 2 £ D r or x\,x 2 £ \ D T . The estimates (E2D and dUSD imply that VX(-,t) extends 

continuously to 8 Dq from either side. I11 the same vein, VX _ 1 (-,t) extends continuously to dD t from 
either side. This follows from standard estimates and Gronwall’s inequality, as we recall below for the 
sake of the reader. One starts by 

X(a,t) = a+ / v(X(oi,t),t) dr. 

Jo 

Using (Ha we obtain 

\X(a,t) -X(P,t)\ <\a- P\+q [ \X(a,r) - X(P,t)\cIt 

Jo 

which yields by Gronwall 


Since 


we get 


and 


\X(a,t)-X(p,t)\ < \a — P\ e qt . 

VX(q, t) = I -f / Vv(X(a, r), r) o VX(q, r) dr, 
Jo 

||VX(-,t )|| 00 < e qt 


\\XX(a,t)-XX(P,t)\\ <q\a-p[< [ e^ +1 ^ dr + [ q\\ VI(a, r) - VI(ft r)|| dr, 

Jo Jo 

provided a, P £ D 0 or a, (5 £ \ D 0 . Thus 


\\VX(a,t)-VX(p,t)\\ < - e ^+^\a-P\y, 

7 


for a, p £ D 0 or a, P £ \ D 0 . 

Since 


—X 1 (x,t) = -v(X %(x,T),t-T) 


the same argument applied to X 1 yields 


\XX-\x,t)-XX-\y,t)\\ < -e^ + ^\x-y\\ 

7 
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for x, y £ D t or x, y £ \ D t . 

For x £ dD t set 

M = M(x,t) = lim VX _1 (y,t) 

DtBy^-x 

and 

N = N(x,t)= lim S/X~ 1 (y,t). 

E d \Dt3y—>x 

Then M(x,t) and N[x,t) are linear mappings from R d into itself that depend continuously on x £ dD t . 
Let Tan(9D 4 ,x) stand for the tangent space to dD t at the point x £ dD t . 

Lemma 8.1. For x £ dD t the linear mappings M and N coincide on Tan {dD t ,x) with the differential 
at x of X~ 1 (x,t), viewed as a differentiable mapping from dD t onto 8 Dq. In particular, M and N map 
Tan (dD t ,x) into Tan(dZ?o,X _1 (x,t)). 

Proof. Let T be a tangent vector to dD t at x, a > 0 and 2 : (—a, a) —> dD t a curve of class C 1 such that 
z(0) = x and z'(0) = T. Let n be the exterior unit normal vector to dD t at x. Let 0 < e < a so small 
that z(9) — rjn £ D t provided \9\ < e and 0 < rj < e. Then 

X~ 1 {z(6) — rjn,t) X~ 1 (z(8),t) uniformly in 0 £ (— e, e). 

Hence 

■^X~ 1 (z(9) - rin,t) X~ 1 {z(9),t) 

as distributions on (—e, e). 

On the other hand, for 6 £ (—e, e) we have 

■^X~ 1 {z(9) - rjn,t) = VI” 1 (z(0) - rjn,t) M(z(9),t) 

pointwise and boundedly. Thus 

±X~\z(9),t) = M(z(9),t) (1^) , 6 £ (—e, e), 

which, for 9 = 0, gives 

DX^IT) = M(T), 

where D.Y -1 is the differential at x of X~ x as a smooth map from dD t into BDq. The argument for N 
is similar. The proof is complete. □ 

Let $0 be a defining function for D 0 . Then $ 0 £ C 1+7 (R d ), D 0 = {x £ R d : <&o(x) < 0}, dD 0 = {x £ 
R d : $o(a ; ) = 0} and Vd>o(x) ^ 0, x £ dD 0 . Set 

ip(x,t) = $o(X~ 1 (x,t)), x £ R d , 

so that D t = {x £ R d : tp(x,t) < 0}, dD t = {x £ : ip(x,t) = 0} and Vtp(y,t) = V$ 0 oVA r ~ 1 (y,t), 

for y (£ dD t . Our next task is computing the jump of Vi p at dD t . 

Let ft and no be the exterior unit normal vectors to dD t and 8Dq respectively, at the points x £ dD t 
and X~ x {x, t) = a £ dD 0 . Given any vector u £ R d , we have 

lim (Vip{y,t),u) = (V$ 0 (a), lim V.Y _1 (?/, t)(u)) = (V$ 0 (a), M(u)). 

D t 3y^x D t 3y^x 
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Set « = An + T, A € R and T G Tan(5D t , x), and 

M(n) = Ano + To, d£l, To £ Tan(<9Do,a). 


Then 


M{u) = XM{n) + M(T) 
and, since A T 0 + M(T ) G Tan(<9D 0 , a), 

(V$ 0 (a), Af(iZ)> = AA(V$ 0 (a),4) 


AAtiq + ATo + M(T) 
AAL|V$ 0 (a)| = A|V$ 0 (a)IM- 


Therefore 

lim V<p(y,t) = A|V$o(aO|n. 

D t By^fX 

Take an orthonormal basis {ti, ...,t^-i} of Tan(SD t ,ir) and an orthonormal basis {rf,..., t^_ ± } of 
Tan(9D 0 ,a) so that 

det(n,ri,.. 1 ) = det(no,r°,.. = 1. 


Call D the differential of X 1 (x, t ) at x as a smooth mapping from dD t into <9-Do- Then the matrix of M 
in the above basis is 

(A 0 ... 0 \ 

Ali 


M = 


Taking determinants 


: Matrix of D 

V Ad-1 

det M = A det D. 


Now detM is the limit of det VX 1 (y,t ) as y G D t tends to x , which turns out to be e*. This is so 
because det VX(a,l) = e -< for a G Dg, which in turn is due to the the well-known fact that (j33] p. 5]) 


dt 


det VX(a, t) 


di vv(X(a,t),t) detVA(a,t) 


det VA(a,f). 


Therefore 


lim 

D t 3y -> 


X<p{y,t) = e‘ 


|V$ 0 (a)| _ 

—;-—— n. 

det D 


Arguing similarly for the exterior side, where det VX 1 (y, t) is 1, we conclude that 


Jim 

R d \Dt3y—>x 


|v^o(q)U 

detD n ' 


Then, clearly, the function 


®(y, t) = e f tp(y, t)X Dt (y) + <p(y, t)X Rd ^{y) 

has no gradient jump at dD t and so is a defining function for D t of class C 1+1 . The material 

derivative of ${y,t) is 

Di =-*(3J,t)X Dt {y) 

and $(y, 0) = $o(j/)- Hence the function determined by (17.21) is of class C 1+1 , as desired. 

Acknowledgements. A.Bertozzi was partially supported by NSF grants CMMF1435709 and DMS- 
0907931 ; J. Garnett by NSF grant DMS 1217239; T. Laurent by NSF DMS 1414396; J. Verdera by 
MTM2013-44699-P (MINECO) and 2014SGR75 (Generalitat de Catalunya). J. Verdera is grateful to the 
Department of Mathematics of UCLA for their hospitality in October 2012. 


27 







References 

[1] Luigi Ambrosio, Edoardo Mainini, and Sylvia Serfaty. Gradient flow of the Chapman-Rubinstein- 
Schatzman model for signed vortices. Annales de I’Institut Henri Poincare (C) Non Linear Analysis, 
28(2):217 - 246, 2011. 

[2] Luigi Ambrosio and Sylvia Serfaty. A gradient flow approach to an evolution problem arising in 
superconductivity. Comm. Pure Appl. Math., 61(11): 1495—1539, 2008. 

[3] N.Baouri, J.Y. Chemin and R.Danchin. Fourier Analysis and Non-Linear Partial Differential Equa¬ 
tions, volume 343 of Studies in Math.. Springer-Verlag, Berlin and Heidelberg, 2011. 

[4] J. Bedrossian, Nancy Rodriguez, and Andrea L. Bertozzi. Local and global well-posedness for aggre¬ 
gation equations and Patlak-Keller-Segel models with degenerate diffusion. Nonlinearity, 24:1683- 
1714, 2011. 

[5] D. Benedetto, E. Caglioti, and M. Pulvirenti. A kinetic equation for granular media. RAIRO Model. 
Math. Anal. Numer., 31 (5):615-641, 1997. 

[6] A. L. Bertozzi and P. Constantin. Global regularity for vortex patches. Comm. Math. Phys., 
152(1):19—28, 1993. 

[7] A. L. Bertozzi, T. B. Laurent, and F. Leger. Aggregation and spreading via the Newtonian potential: 
the dynamics of patch solutions. Math. Models and Methods in the Applied Sciences, 22(Supp. 01), 
2012 . 

[8] A.L. Bertozzi and T. Laurent. Finite-time blow-up of solutions of an aggregation equation in R n . 
Comm. Math. Phys., 274:717-735, 2007. 

[9] Andrea L. Bertozzi and Jeremy Brandman. Finite-time blow-up of L°°-weak solutions of an aggre¬ 
gation equation. Commun. Math. Sci., 8(l):45-65, 2010. 

[10] Andrea L. Bertozzi, Jose A. Carrillo, and Thomas Laurent. Blow-up in multidimensional aggregation 
equations with mildly singular interaction kernels. Nonlinearity, 22(3):683-710, 2009. 

[11] Andrea L. Bertozzi, John B. Garnett, and Thomas Laurent. Characterization of radially symmet¬ 
ric finite time blowup in multidimensional aggregation equations. SIAM Journal on Mathematical 
Analysis, 44(2):651-681, 2012. 

[12] Andrea L. Bertozzi, Thomas Laurent, and Jesus Rosado. L p theory for the multidimensional aggre¬ 
gation equation. Comm. Pure Appl. Math., 64(l):45-83, 2011. 

[13] Andrea L. Bertozzi and Dejan Slepcev. Existence and uniqueness of solutions to an aggregation 
equation with degenerate diffusion. Commun. Pure Appl. Anal., 9(6):1617-1637, 2010. 

[14] Piotr Biler and Wojbor A. Woyczynski. Global and exploding solutions for nonlocal quadratic 
evolution problems. SIAM J. Appl. Math., 59(3):845-869 (electronic), 1999. 


28 


[15] Adrien Blanchet, Jose A. Carrillo, and Nader Masmoudi. Infinite time aggregation for the critical 
Patlak-Keller-Segel model in R 2 . Comm. Pure Appl. Math., 61(10):1449-1481, 2008. 

[16] Adrien Blanchet, Jean Dolbeault, and Benoit Perthame. Two-dimensional Keller-Segel model: op¬ 
timal critical mass and qualitative properties of the solutions. Electron. J. Differential Equations, 
pages No. 44, 32 pp. (electronic), 2006. 

[17] M. Bodnar and J. J. L. Velazquez. An integro-differential equation arising as a limit of individual 
cell-based models. J. Differential Equations, 222(2):341 380, 2006. 

[18] Silvia Boi, Vincenzo Capasso, and Daniela Morale. Modeling the aggregative behavior of ants of 
the species polyergus rufescens. Nonlinear Anal. Real World Appl, 1(1):163-176, 2000. Spatial 
heterogeneity in ecological models (Alcala de Henares, 1998). 

[19] Michael P. Brenner, Peter Constantin, Leo P. Kadanoff, Alain Schenkel, and Shankar C. Venkatara- 
mani. Diffusion, attraction and collapse. Nonlinearity, 12(4): 1071—1098, 1999. 

[20] Martin Burger, Vincenzo Capasso, and Daniela Morale. On an aggregation model with long and 
short range interactions. Nonlinear Anal. Real World Appl., 8(3):939-958, 2007. 

[21] Martin Burger and Marco Di Francesco. Large time behavior of nonlocal aggregation models with 
nonlinear diffusion. Netw. Heterog. Media, 3(4):749-785, 2008. 

[22] J. A. Carrillo, M. DiFrancesco, A. Figalli, T. Laurent, and D. Slepcev. Global-in-time weak measure 
solutions and finite-time aggregation for nonlocal interaction equations. Duke Math. J., 156(2):229- 
271, 2011. 

[23] Jose A. Carrillo, Robert J. McCann, and Cedric Villani. Kinetic equilibration rates for granular 
media and related equations: entropy dissipation and mass transportation estimates. Rev. Mat. 
Iberoamericana, 19(3):971 1018, 2003. 

[24] Jose A. Carrillo, Robert J. McCann, and Cedric Villani. Contractions in the 2-Wasserstein length 
space and thermalization of granular media. Arch. Ration. Mech. Anal., 179(2) :217-263, 2006. 

[25] Jose A. Carrillo and Jesus Rosado. Uniqueness of bounded solutions to aggregation equations by 
optimal transport methods. In European Congress of Mathematics, pages 3-16. Eur. Math. Soc., 
Zurich, 2010. 

[26] Angel Castro, Diego Cordoba, Charles Ferfferman and Francisco Gancedo. Breakdown of Smoothness 
for the Muskat Problem. Arch. Ration. Mech. Anal., 208(3):805-909, 2013. 

[27] Angel Castro, Diego Cordoba, Charles Ferfferman, Francisco Gancedo and Javier Gomez-Serrano. 
Finite time singularities for the free boundary incompressible Euler equations. Annals of Math., 
178(3):1061 1134,2013. 

[28] Jean-Yves Chemin. Persistance de structures geometriques dans les fluides incompressibles bidimen- 
sionnels. Ann. Sci. Ecole Norm. Sup. (4), 26(4):517-542, 1993. 


29 


[29] Diego Cordoba and Francisco Gancedo. Absence of Squirt Singularities for the Multi-Phase Muskat 
Problem. Comm. Math. Phys., 299(2/561-575, 2010. 

[30] Jean Dolbeault and Benoit Perthame. Optimal critical mass in the two-dimensional Keller-Segel 
model in R 2 . C. R. Math. Acad. Sci. Paris , 339(9):611-616, 2004. 

[31] Hongjie Dong. On similarity solutions to the multidimensional aggregation equation. SIAM J. Math. 
Anal, 43(4):1995-2008, 2012. http://arxiv.org/abs/1102.0177. 

[32] Qiang Du and Ping Zhang. Existence of weak solutions to some vortex density models. SIAM J. 
Math. Anal., 34(6/1279-1299 (electronic), 2003. 

[33] Weinan E. Dynamics of vortex liquids in Ginzburg-Landau theories with applications to supercon¬ 
ductivity. Physical Review B, 50(2): 1126—1135, 1994. 

[34] Veysel Gazi and Kevin M. Passino. Stability analysis of swarms. IEEE Trans. Automat. Control, 
48(4):692-697, 2003. 

[35] Darryl D. Holm and Vakhtang Putkaradze. Formation of clumps and patches in self-aggregation of 
finite-size particles. Phys. D, 220(2):183—196, 2006. 

[36] Yanghong Huang and Andrea L. Bertozzi. Self-similar blowup solutions to an aggregation equation 
in R". SIAM J. Appl. Math., 70(7):2582-2603, 2010. 

[37] E. F. Keller and L. A. Segel. Initiation of slime mold aggregation viewed as an instability. J. Theor. 
Biol., 26:399-415, 1970. 

[38] Thomas Laurent. Local and global existence for an aggregation equation. Comm. Partial Differential 
Equations, 32(10-12/1941 1964, 2007. 

[39] Dong Li and Jose Rodrigo. Finite-time singularities of an aggregation equation in R™ with fractional 
dissipation. Comm. Math. Phys., 287(2):687-703, 2009. 

[40] Dong Li and Jose L. Rodrigo. Refined blowup criteria and nonsymmetric blowup of an aggregation 
equation. Adv. Math., 220(6/1717-1738, 2009. 

[41] Dong Li and Xiaoyi Zhang. On a nonlocal aggregation model with nonlinear diffusion. Discrete 
Contin. Dyn. Syst., 27(1/301-323, 2010. 

[42] Hailiang Li and Giuseppe Toscani. Long-time asymptotics of kinetic models of granular flows. Arch. 
Ration. Mech. Anal., 172(3/407-428, 2004. 

[43] Andrew J. Majda and Andrea L. Bertozzi. Vorticity and incompressible flow, volume 27 of Cambridge 
Texts in Applied Mathematics. Cambridge University Press, Cambridge, 2002. 

[44] Nader Masmoudi and Ping Zhang. Global solutions to vortex density equations arising from sup- 
conductivity. Ann. Inst. H. Poincare Anal. Non Lineaire, 22(4/441-458, 2005. 


30 



[45] Joan Mateu, Joan Orobitg, and Joan Verdera. Extra cancellation of even Calderon-Zygmund 
operators and quasiconformal mappings. Journal de Mathematiques Pures et Appliques, 91(4):402- 
431, 2009. 

[46] Joan Mateu, Joan Orobitg, and Joan Verdera. Estimates for the maximal singular integral in 
terms of the singular integral : the case of even kernels. Annals of Math., 174(3):1429-T483, 2011. 

[47] Alexander Mogilner, Leah Edelstein-Keshet, L. Bent, and A. Spiros. Mutual interactions, potentials, 
and individual distance in a social aggregation. J. Math. Biol., 47(4):353-389, 2003. 

[48] Alexander Mogilner and Leah Edelstein-Keshet. A non-local model for a swarm. J. Math. Biol., 
38(6) :534—570, 1999. 

[49] Daniela Morale, Vincenzo Capasso, and Karl Oelschlager. An interacting particle system modelling 
aggregation behavior: from individuals to populations. J. Math. Biol., 50(1):49—66, 2005. 

[50] Juan Nieto, Frederic Poupaud, and Juan Soler. High-field limit for the Vlasov-Poisson-Fokker-Planck 
system. Arch. Ration. Mech. Anal., 158(1):29—59, 2001. 

[51] Frederic Poupaud. Diagonal defect measures, adhesion dynamics and Euler equation. Methods Appl. 
Anal., 9(4):533—561, 2002. 

[52] Silvya Serfaty and Juan Luis Vazquez. A mean field equation as limit of non-linear diffusions with 
fractional Laplacian operators. Calc. Var., 49:1091-1120, 2014. 

[53] Chad M. Topaz and Andrea L. Bertozzi. Swarming patterns in a two-dimensional kinematic model 
for biological groups. SIAM J. Appl. Math., 65(1):152-174 (electronic), 2004. 

[54] Chad M. Topaz, Andrea L. Bertozzi, and Mark A. Lewis. A nonlocal continuum model for biological 
aggregation. Bull. Math. Biol, 68(7):1601 T623, 2006. 

[55] Giuseppe Toscani. One-dimensional kinetic models of granular flows. M2AN Math. Model. Numer. 
Anal, 34(6):1277-1291, 2000. 

[56] V. I. Yudovich. Non-stationary flow of an ideal incompressible liquid. USSR Computational Mathe¬ 
matics and Mathematical Physics, 3(6):1407 1456, 1963. 


31 


Andrea L. Bertozzi 
Department of Mathematics 
University of California at Los Angeles 
Los Angeles, CA 90095 

E-mail: bertozzi@math.ucla.edu 

John B. Garnett 
Department of Mathematics 
University of California at Los Angeles 
Los Angeles, CA 90095 

E-mail: jbg@math.ucla.edu 

Thomas B. Laurent 
Department of Mathematics 
Loyola Marymount University 
Los Angeles, CA 90045 

E-mail: tlaurent@lmu.edu 

Joan Verdera 

Departament de Matematiques 
Universitat Autonoma de Barcelona 
08193 Bellaterra, Barcelona, Catalonia 

E-mail: j vm@mat. uab. cat 


32 


